Weak-Lensing by Large-Scale Structure and the 
Polarization Properties of Distant Radio- Sources 



Gabriela C. Surpi and Diego D. Harari 
Departamento de Fisica, Facultad de Ciencias Exactas y Naturales 
Universidad de Buenos Aires 
Ciudad Universitaria - Pab. 1, 1428 Buenos Aires, Argentina 

ABSTRACT 

We estimate the effects of weak lensing by large-scale density inhomogeneities 
and long-wavelength gravitational waves upon the polarization properties of 
electromagnetic radiation as it propagates from cosmologically distant sources. 
Scalar (density) fluctuations do not rotate neither the plane of polarization of the 
electromagnetic radiation nor the source image. They produce, however, an appreciable 
shear, which distorts the image shape, leading to an apparent rotation of the image 
orientation relative to its plane of polarization. In sources with large ellipticity 
the apparent rotation is rather small, of the order (in radians) of the dimensionless 
shear. The effect is larger at smaller source eccentricity. A shear of 1% can induce 
apparent rotations of around 5 degrees in radio sources with the smallest eccentricity 
among those with a significant degree of integrated linear polarization. We discuss 
the possibility that weak lensing by shear with rms value around or below 5% may 
be the cause for the dispersion in the direction of integrated linear polarization of 
cosmologically distant radio sources away from the perpendicular to their major axis, 
as expected from models for their magnetic fields. A rms shear larger than 5% would 
be incompatible with the observed correlation between polarization properties and 
source orientation in distant radio galaxies and quasars. Gravity waves do rotate 
both the plane of polarization as well as the source image. Their weak lensing effects, 
however, are negligible. 



Subject headings: cosmology: theory — gravitational lensing — large-scale structure 
of universe — polarization 



Introduction 



Polarization properties of extragalactic radio sources are primarily a source of information 
about magnetic fields in the source, in the intergalactic medium, and in the Milky Way ( ^aikia| 



fe Salter 1988, Kronberg 1994). The orientation of the plane of linear polarization at the source, 
derived from multi-frequency polarization measurements after correction for Faraday rotation 
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along intervening magnetic fields, should point perpendicularly to the magnetic field lines, the 
emission being synchrotron radiation. 

In the case of extended radio jets from quasars, the orientation of the intrinsic linear 
polarization within a well defined jet is typically perpendicular to the local direction of each jet 



segment, and if it does change it is usually a 90° flip (see, e.g., Bridle, Per ley &: Henriksen 1986 ; 
and [Killen, Bicknell, fc Ekers 1986| ). Computational MHD modeling of radio jets improves the 
confidence in the predicted jet-intrinsic magnetic field configurations supporting these observations 



( Clarke, Norman, fc Burns 198£ 



Kronberg et al. (1991) have proposed and applied a technique involving the angular 
relationship between the intrinsic polarization vectors and the morphological structure of extended 
radio jets to probe the gravitational distortion of these sources by foreground lensing galaxies. 
Since the polarization vector is not rotated by the gravitational field of an ordinary gravitational 
lens, the bending of a jet away from its intrinsic projected shape caused by the lensing effects of 
intervening masses yields a detectable local "alignment breaking" between the polarization and 
the orientation of each jet segment. Application of this idea to the analysis of the radio and deep 
optical images of the strongly polarized jet from 3C 9 provided the first estimate of the mass and 
mass distributions of two intervening galaxies acting as lenses using this technique (Kronberg et 
al. 1991 ; and Kronberg, Dyer, &: Roser 1996| ). 



The notion that the plane of polarization is not rotated by the gravitational field of an 
ordinary gravitational lens, which underlies the method of Kronberg et al. (1991), was addressed 
by Dyer &; Shaver (1992). Using symmetry considerations they argued that the polarization 
direction of photons in a beam under the influence of a gravitational lens remains unaffected 
except for lenses in very rapid rotation, an astrophysically unlikely situation. Similar conclusions 
were also reached by Faraoni (1993), with a different approach. The fact that the gravitational 
field of rotating bodies may rotate the plane of polarization of electromagnetic waves has been 
addressed among others by Strotskii (1957), Balazs (1958), Plebanski (1960), Mashhoon (1973), 
Pinneault & Roeder (1977) and Su & Mallet (1980). 

In the present work we analyse the effect of weak gravitational lensing by large scale structure 
in the Universe upon the polarization properties of electromagnetic radiation from distant sources. 
Through an extension of the method of "alignment breaking" , introduced by Kronberg et al. 
(1991), we suggest that polarization properties of distant radio sources may provide information 
on the fluctuations in the gravitational potential due to large scale density inhomogeneities. 

When light from a distant source propagates across density inhomogeneities, the apparent 
shape and size of the source is distorted by tidal deflections. The effect is known as weak lensing, 
in contrast to the case when the bending of light rays is so strong that multiple image formation 
is possible. Weak lensing effects have already proved to be a powerful tool to map the mass 



distribution in galactic halos (Tyson et al. 1984) and clusters of galaxies (Tyson, Valdes, & Wenk 



1990). They are also a potential probe of density fluctuations on the largest scales ( [BlandfoTd 
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et al. 1991, [Miralda-Escude 1991| , ^ aiser 199^ ), particularly through the correlation in galactic 



ellipticities induced by the shear they produce. It has been estimated that shear caused by 
large-scale density fluctuations in the linear regime is of the order of 1%, depending on the 
cosmological model, the normalization of the power spectrum, and the redshift distribution of 
sources (Blandford et al. 1991, [Miralda-Escude 1991 , Kaiser 1992). Theoretical predictions on 



smaller scales (below ^^20 arcmin) by Jain & Seljak (1996), show that density inhomogeneities 
in the non-linear regime induce a larger cosmic shear, of the order of a few percent. Recently, 
Schneider et al. (1997) presented the first evidence for the detection of a robust shear signal of 
about 3% on scales of 1 arcmin, in agreement with the above theoretical expectations, from the 
analysis of data in a field containing the radio source PKS 1508-05. Long- wavelength gravitational 
waves (tensor metric perturbations) are also a potential source of weak lensing of distant sources. 

Radio sources with a significant degree of linear polarization, are also usually elongated along 
one direction. The angle x ~ between the intrinsic plane of polarization and the orientation of 
the structural axis has been measured for a large number of sources, for which redshift information 
is also available. The existence of a significant correlation between the direction of integrated 
linear polarization and the source axis orientation has been noticed since the earliest investigations 



( Haves &: Conway 1975 , Clarke et al. 1980 ). For high redshift galaxies, a well defined peak is 
found around X ~ = 90° > suggesting a magnetic field parallel to the source axis. For low redshift 
sources there is weaker correlation, with a narrow peak at x — ^ = 0°, suggesting a magnetic 
field perpendicular to the source axis, and a broader peak at x ~ V' = 90°- It is believed that the 
X — V' = 90° peak corresponds to high-luminosity sources, while the x ~ V' = 0° peak is ascribed to 
low-luminosity objects. 

The angle between the plane of integrated linear polarization and the major axis of distant 
radio sources is potentially a useful tool to search for unusual effects upon light propagation over 
cosmological distances, additional to the standard Faraday rotation. There have been, for instance, 
searches for correlations between the polarization properties of extragalactic radio sources and 
their location, which would evidence some kind of "cosmological birefringence". While first 



attempts indicated that the residual polarization rotation fitted a dipole rule ( pirch 1982| , [Kendall 



[fc Young 1984| , Bietenholz fc Kronberg 1984| ), a later analysis of a larger sample by Bietenholz 



(1986) rejected that evidence. Recently, Nodland & Ralston (1997) have claimed to find evidence 
in the data for a redshift-dependent dipole anisotropy in the rotation of the plane of polarization. 
A more conventional interpretation and statistical analysis of the same data ( Carroll fc Field] 



1997, Eisenstein fc Bunn 1997] ) rejects any statistically significant positive signal, which would 



also be in conflict with other potentially more accurate probes of such cosmological rotation of 
polarization, through high resolution polarization and intensity data of some specific radio sources 
([Leahy 1997[ , [Wardle et al. WdTj ). 



In this paper, we estimate the effect of weak lensing by large scale structure, in the form 
of density fluctuations or gravitational waves, upon the observed angle between the plane of 
polarization of the electromagnetic emission by distant radio sources and their image orientation. 
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We find it is in principle possible that the scattered cases of distant radio sources which appear 
to be polarized away from the direction either perpendicular or parallel to the source axis, 
if Faraday rotation were the only propagation effect, be actually a consequence of additional 
wavelength-independent rotation of the polarization as their light travels to us through an 
inhomogeneous universe. 

In section 2 we review the formalism appropriate to describe image distortions caused by 
scalar and tensor weak lensing. In section 3 we work in the geometric optics approximation and 
derive the rotation of the plane of polarization as light propagates in the presence of metric 
fluctuations. In section 4 we estimate the effect of large scale structure upon the angle x ~ ^ 
between the plane of polarization and the orientation of a source image, and show that a shear of 
1% can change the value of x — by 5°. In section 5, we discuss the possibility that scalar shear 
with a rms value around or below 5%, which would cause apparent rotations with a mean square 
dispersion of 20°, be the cause of the scattered cases of cosmologically distant radio sources with 
a direction of integrated linear polarization which is not perpendicular to their major axis. In 
section 6 we summarize our conclusions. 



2. Image distortion by random metric perturbations 



In this section we review the formalism appropriate to describe weak lensing effects upon 
light propagation along cosmological distances through a universe with large scale density 
inhomogeneities and long wavelength gravitational waves. Part of the material in this section is 
not new (see for instance Kaiser 1992 , ^eljak 1994 , Kaiser fc Jaffe 1997 and Kaiser 1996 ). We 
include it for completeness, to define our method and conventions, and because some of our results 
for the gravitational distortions are at variance with those of Seljak (1994) and Kaiser & Jaffe 
(1997). The differences do not conflict with the main conclusion that weak lensing effects by large 
scale density inhomogeneities accumulate over distances longer than the wavelength of the metric 
perturbation as (Z)/A)^/^ ( peljak 1994 ), while weak lensing distortions by gravitational waves do 
not significantly accumulate ( Kaiser fc Jaffe 1997] ). Even if the effect of gravitational waves is 
small, it is relevant to precisely determine the focusing, shear and rotation they induce if we want 
to determine the rotation of the plane of polarization relative to the morphological structure of 
the source. 

We restrict our attention to effects that are linear in the metric perturbations. For 
simplicity, we shall consider scalar and tensor metric fluctuations around a Minkowski background. 
The generalization to a Robertson- Walker background is straightforward ( Kaiser 1996| ). The 
restriction to a Minkowski background oversimplifies the distance-dependence of the effects under 
consideration, but preserves their qualitative features as well as the order of magnitude estimates. 
Vector perturbations could also be easily incorporated into this formalism. 



Consider a Minkowski background metric r]^y with signature (-1,1,1,1), and linear 
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perturbations Greek indices nv,a,(3,... = 0,1,2,3 denote space-time coordinates. Latin 
indices = 1,2,3 will denote spatial coordinates, but latin indices a,b,... = 1,2 will be 

reserved for the components transverse to the unperturbed photons paths, that we shall take as 
approximately parallel to the x^— axis. Boldface symbols will denote spatial vectors, i.e. m, k,... 

The linearized geodesic equations read, 

^ = - ^ • (1) 

We denote the photons 4-momentum with = dx^/dz, where z is an affine parameter (not 
necessarily coincident with the x'^-coordinate) . 

We consider photon paths sufficiently close to the x^— axis, pointing towards an observer 
located at the origin of coordinates, and work not only up to linear order in the metric 
perturbations, but also up to first order in the photon departures away from the polar axis, which 
we parametrize as 

x^iz) = 9''z + 0{h,e^) (a = 1,2) . (2) 

We choose an affine parameter such that x^ = z + O{h,0'^), and thus t = + {zs — z) + 0{h, 0'^), 
with tg the time at which the photon is emitted at the source, located at an affine distance Zs in 
the x^-direction. We assume a gauge such that hoj = 0. To lowest order p^^ = {—1,9^,9'^, 1), and 
the geodesic equation for the transverse departure of the photon path away from the x^— axis can 
be written as: ^ 

D''{^,t) + F\{^,t)9^ + 0{9^) (3) 



dz^ 

where 



Da — ^{hoo + hss)^a + ha3,0 — ha3,3 ^^-j 
Fab = habfi — hab,3 — {h3a,b — h^h^a) ■ 



Taylor-expanding Da{:x.,t) to first order in x" we can finally write 



D^iz) + M''h{z)9' + 0(9') (5) 



dz^ 

where we defined 

Mab{z) = zDa,b{z) + Fab{z) . (6) 

The quantities in the r.h.s of equation (||) are evaluated over the x^-axis and along the unperturbed 
photon path. For instance Da{z) = Da{x°- = z^t = + z^ — z). 

In order to determine the distortion produced by metric fluctuations upon light propagation, 
we consider two neighboring photon paths that arrive to the observer with angular separation 0" 
coming from two different points at the source separated by a coordinate distance /S.x"'[zs) = 9'^Zs. 
The particular solution to the geodesic deviation equation with the appropriate boundary 
conditions is ^ 

Ax"(z) = rz+ r dz' r dz"M\{z")9^ . (7) 
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A source with a shape described in these coordinates by Ax^ = O'^Zg is seen by an observer at the 
origin as if the image were at Ax'^ = O'^Zg in the source sphere. The mapping 9"" — > 0" is given by 

e'' = {5\ + rb)0' (8) 

where 

i^ab = - r dz^^Mahiz) (9) 

Jo Zs 

is the distortion matrix, which describes the effect of hnear metric perturbations upon the observed 
shape of an image. 

We point out now what is the source of variance between our results and those by Seljak 
(1994) and Kaiser & Jaffe (1997). Their expressions for the geodesic deviation and ipab coincide 
with eqs. (0) and (P) with Mat = zDa^b, rather than the complete expression in eq. (P). The 
difference is the term Fat, which arises as a consequence that the photon paths must converge into 
the origin, and thus the unperturbed paths are not strictly parallel. The relevant observational 
situation is not exactly described by the distortion suffered by a bundle of initially parallel photon 
paths, but rather by a bundle that starting at the source focalize at the observer location. In 
other words, tidal deflections are not only due to the gradients of the gravitational potential in 
the directions transverse to the central ray in the bundle, but there is also a contribution of the 
gradient in the longitudinal direction, given that the rays in the bundle are not parallel. The 
longitudinal gradients, as we shall see, are significant in the tensor case only, which at any rate is 
not likely to have a large impact upon measurements. The precise result for the focusing, shear 
and rotation induced by gravity waves, however, must take them into account. 

There are some subtleties in the definition of the distortion matrix to linear order in the 
perturbations that we now discuss. As defined by eq. (P), the distortion matrix relates coordinate 
(rather than proper) shapes. It is thus convenient to separate from jpab the local distortion effects, 
due to the different proper length that a same coordinate segment has at the observer and source 
locations. We thus define ^pab such that 

i^ab = i'ab + ^Ahab ■ (10) 

Here Ahab = hab{zs,te) — hab{z = 0,t = to) is the difference in the metric perturbations between 
the events of emission of the light ray at time te from the source location, and of its observation 
at the origin at time to- ipab is the appropriate matrix to describe the image distortion in terms of 
proper length mappings: 

gab{0, Qe'^e' zl = [gab{Zs,t,) + i^ab] e^O' zl (11) 

with gab = rjab + Kb- 

Although this is just a conventional matter, we also wish to emphasize that our definition of 
focusing describes its dependence with affine, rather than coordinate, distance. Since they differ 
by terms linear in the perturbations, this matter of convention needs to be specified. 
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The more general distortion matrix ipab can be decomposed in terms of a diagonal matrix, 
which represents focusing or convergence, a symmetric traceless matrix 7^^, which describes shear, 
and an antisymmetric matrix, that corresponds to a rigid rotation of the image: 

Ipab = Kdab + ^ab - ^(-ab ■ (12) 

Here eab is the Levi-Civita tensor in two-dimensions. Shear can be described in terms of its 
intensity 7 and the direction of one of its principal axis (j) (which is only defined modulo vr) as 

/ cos2(/) sin 20 \ 

7afe = 7 . „ , „ , ■ (13) 

\ sm 20 — cos 20 / 

Notice that a rotation of the direction of the shear by 7r/2 amounts to just a change in the sign of 
7. Two independent directions in which to decompose an arbitrary shear form an angle 7r/4. 

We now specialize the expressions above to the case of scalar and tensor metric perturbations 
respectively. We start working with a single planar metric fluctuation, and later compute the rms 
effect of a full stochastic background by superposition. 



2.1. Scalar perturbations 

We consider one Fourier mode of static scalar perturbations 

V(x) = h{\^)e'''-^5^p = h{\^)e'^^^\-'^^'-^5^p (14) 

appropriate for instance to describe adiabatic density perturbations in a matter-dominated, 
spatially-flat Robertson- Walker background, in the longitudinal gauge. We have defined ^ = m • n, 
where /i is the cosine of the angle between the wave-vector k = km, and the direction n in which 
light propagates, which is approximately along the x^— axis, incoming towards the origin from 
positive values of . Specialization of the eqs. in the previous section to this case lead to 

ipab = h{\<i)[a{ii,Zs)mamh + f5{ii,Zs)li5ab] ■ (15) 

Here we defined 

a(^,z.) = e dz'-^e-^^^^ ^^g^ 



PipL.Zs) = -ik J^' dz^^^e-^''^''' . 

We write the projection of the wave- vector direction onto the plane orthogonal to z as 

m"' = y/l — /x^(cos (f, sin ip), and from eqs. ([l2|) and (15) we find that the focusing, shear intensity 

and direction, and rotation induced by a single Fourier mode of scalar perturbations are 

K = i/i(k) [(1 - iJ?)a{iJi, Zs) + 2fif3{n, Zs)] 

(17) 



i/i(k) (1 - iJ?)a{iJL,Zs) 



w = 
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Scalar perturbations induce no rotation. They produce shear in the direction of the perturbation 
wave-vector, projected into the source plane. 

As pointed out by Kaiser & Jaffe (1997), the reason why weak lensing effects by scalar 
perturbations accumulate over distances longer than the perturbation wavelength is evident 
from eqs. ( p^ ) and (|l^). Resonant modes, those for which the photons remain in phase with 
the perturbation, are those with ^ ~ 0, which produce a significant deflection. Notice that 
a{iJ,) {kzsf' and /3(/i) — > ikzs/l as — > 0. Our results differ from those by Seljak (1994) 
and Kaiser & Jaffe (1997) by the effect of the longitudinal modes which affect just the focusing 
K by the term proportional to fj,(3{fi, Zg). Since this term does not contribute in the resonant 
situation (^ = 0), it is sub-dominant, and negligible, after superposition of all Fourier modes of 
the stochastic background. 

As discussed in the previous subsection, the local effects can be subtracted from the 
distortion matrix describing proper distortions through ipab = ipab — ^^^afe- In this case 
Ahab = 5ab{e~^^^^" — l)i and use of the relation 2/u/3 — /i^a = (e~*'^^^« — 1) allows to identify the 
proper distortion matrix elements as 

K = i/i(k)a(/x,Zs) ^^g^ 
7=7- 

Consider now a stochastic background of scalar fluctuations characterized in momentum-space 
by the correlations < /i(k)/i*(k') >= (27r)^(5(k — ]<i!)P{k) where P{k) is the power spectrum. The 
mean square expectation value of the shear induced is 

In the limit kzg ^ 1 the /x-integral equals {kzsY + 0{{kzs)'^\. Shear induced by scalar 
perturbations of amplitude h accumulates over distances D longer than the wavelength A as 
h X (D/A)^/^ (Seljak 1994). In the case of a scale invariant power spectrum P{k), such as that 
predicted by inflationary models in a CDM scenario, with amplitude, normalized at small k 
according to the COBE-DMR measurement of the cosmic microwave background anisotropy, the 
rms shear is of order < 7^ >^/^ss 0.02. The rms convergence is of the same order of magnitude. 



2.2. Tensor perturbations 

Now we study image distortions produced by gravity waves. We work tensor perturbations in 
the transverse traceless gauge, where h^pk"' = 0, /lo/^ = 0, and = 0. We consider gravitational 
waves with wave vector k^ = A;(l,m), with m* = (-\/l — cos 93, ^1 — sin 93, — ^) . As in the 
previous section, /i = m • n is the cosine of the angle between the perturbation wave-vector 
and the photon propagation direction. We denote the two independent polarization modes of 
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the gravity waves with the symbols +, x. The normal modes of the metric perturbations read 
hj{yL,t) = hij{k)e^^^-''-''^'^ with 



/iii(k) = /i4_(k)(/x^ cos^ if — sin^ ip) — hx (k)/i sin 2ip 

^22 (k) = /i+(k)(/u^ sin^ (/? — cos^ ip) + hx (k)// sin 2ip 

/i33(k) =/l+(k)(l-^2) 

/ii2(k) = /i+(k)i(l + fi^)s'm2ip + hxO^jfJ- cos 2(p 

/ii3(k) = /i_|_(k)^-\/l — ij? cosip — /ix(k)-\/l — /i^sin(/3 

/i23(k) = /i+(k)^v/l - ijflsimp + /ix(k)A/l - ijflcosif 

The distortion matrix now reads, from eq. @ 



(20) 



■006= a(;U-l,2;s) 



^hsslTlamb - (1 - fJ-)h3amb - 
(1 - ^)/la6 + (/l3a"lfe - /l3fe"la)] 



(21) 



where a and /3 are exactly the same functions as in eq. (16), now the argument being /i — 1 
instead of fi. We have discarded an irrelevant overall phase factor. The focusing, shear intensity 
and direction, and rotation induced by each normal mode of a gravity wave background, for each 
independent polarization, thus turn out to be 



K =1(1- ^')/i+(k) [(1 - M)'a(/" - 1, zs) + 2(1 - - 1, zs)] 

7+ = i/i+(k) [(1 - _ ^)2c,(^ _ 1, - 2(1 + /i2)(l - ^)l3{is - 1, Zs 

(/)+ =ip 

7x = Ihx (k) [(1 - _ ^)a(/x - 1, Zs) + 2^{1 - ^,)P{^l - 1, z,)] 

UJ =-\{i- i?)hx{\^)[{i- ^J)a{|Ji-l,zs) + 2p{^l-l,zs)] . 



(22) 



Gravity waves have quite different weak lensing effects than scalar metric perturbations. The 
fact that waves propagate imply that the resonant modes, those that may eventually add- up over 
distances longer than the wavelength, are now those with ^ ~ 1, which do not produce neither 
focusing, shear, nor rotation. Eqs. (^) coincide with the equivalent expressions of Kaiser &; Jaffe 
(1997), except for the terms proportional to /3(^ — 1, Zs) 



The + polarization of the gravity waves produces focusing, shear with a principal axis in 
the direction of the wave- vector projection into the source plane {ip), and no rotation. The x 
polarization produces no focusing, a shear comparable in strength to that of the + polarization 
but with principal axis at 45°, and induces rigid rotation of the image. In other words, the 
+ polarization produces similar effects as scalar shear, while the x polarization induces truly 
pseudo-scalar shear (the change from a right-handed coordinate system to a left-handed one 
changes the sign oi hx but not of /i_|_). 

We now use the relation 2(1 — ^)/3 + (1 — //)^a = 1 — express the rotation uj as 

cu = ^(l + ^)A/ix(k) (23) 
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where A/ix(k) = h,,{k){e'''^'^-''>' - 1) is the variation in the amphtude of the x component of the 
gravitational wave between the events of emission of the hght ray at the source and its arrival to 
the observer at the origin. Clearly, there is no possible cumulative effect. The rotation can never 
exceed the dimensionless wave amplitude. 

To separate the local distortion effects from those produced during light propagation we 
evaluate ipab = i'ab — ^^hab, in this case the proper distortion matrix elements read: 

K =0 

7(+,x) = ^(1 - ^)^a(^ - l,Zs)/i(+,x) (24) 

UJ = UJ 

Expressed in terms of proper areas and as a function of affine distance to the source, there is no 
focusing to linear order in the gravity wave amplitude, in agreement with previous conclusions by 
Zipoy & Bertotti (1968). 

Consider now a stochastic background of gravitational radiation, built as a superposition 
of plane sinusoidal waves. Contrary to the scalar case, they induce distortions just of the 
order of their dimensionless amplitude, but not larger (Kaiser & Jaffe 1997). Consider, for 
instance, an stochastic background such as that predicted by inflationary cosmological models 



( Abbott Sz Wise 1984 ), (rms amplitude proportional to the wavelength), with the largest possible 
amplitude compatible with the observed CMB anisotropy. The rms amplitude is of the order of 
h{k) ~ 10~^Ho/k. The dominant effect is now that of the wavelengths comparable to the distance 
to the source. Distortions of order 7 ~ 10"^, a factor 10~^ smaller than the effect expected to 
be induced by scalar perturbations, are unlikely to be detected. At any rate, it is interesting 
to observe that there are specific footprints of weak lensing by tensor perturbations, due to its 
pseudo-scalar contribution, which may eventually serve to distinguish it from that produced by 
density inhomogeneities. 



3. Rotation of the plane of polarization by gravity waves 

One specific footprint of weak lensing by gravity waves is rotation, which is not induced by 
scalar perturbations. Rotation of an image may be difficult to ascertain, without a knowledge of 
the intrinsic source orientation. The polarization properties of the electromagnetic emission of 
the source may, in some cases, bear some correlation with the intrinsic source orientation. The 
angle between some structural direction in a source and the direction of the plane of polarization 
is thus potentially a useful tool to search for eventual image rotation by weak lensing. Since 
the plane of polarization may also rotate as electromagnetic radiation propagates across density 
inhomogeneities or gravitational waves, we should also estimate this effect. To do so, in this section 
we work within the geometric optics approximation of the solution to Maxwell equations in curved 
space-time ( Misner, Thorne, &: Wheeler 1973| ) , which imply that the direction of linear polarization 



is parallel-transported along the photon path. It is intuitive that scalar perturbations, which have 
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no handedness, do not rotate the plane of polarization. In the case of gravity waves, however, we 
may expect them to induce a "gravitational Faraday effect" ( perlick fc Hasse 1995 , Cooperstock fc 



Faraoni 1993| ). Image rotation being a consequence of the geodesic deviation between neighboring 



rays while polarization rotation arises from parallel transport, we may in principle expect them to 
be different, and thus potentially provide a test to detect absolute rotations. 

Typical wavelenghs of electromagnetic radiation are extremely short compared to the radius 
of curvature of space-time and to the typical length over which the amplitude, polarization, and 
wavelength of the electromagnetic waves may vary. Thus, we may confidently apply the geometric 
optics approximation in our analysis, and electromagnetic radiation may be regarded locally as 
plane waves that propagate through a space-time of negligible curvature. In the Lorentz gauge the 
electromagnetic vector potential may be written as 

A^' = ^{af'e'^} (25) 

where is a vector amplitude and is a real phase that rapidly changes, with the typical 
frequency of the electromagnetic radiation, 9 = uop^^x^ + const., with uj the wave frequency and 
the dimensionless wave-vector. Expanding the solution of Maxwell's equations in powers of the 
electromagnetic wavelength divided by the length-scale over which the gravitational potentials 
vary significantly (the gravitational wave wavelength in our case) one finds to lowest order that 
the photon paths are null geodesic, and to the next order that the vector amplitude can be 
written as = af^ where a is a scalar amplitude and is the polarization vector (normalized 
^-s f^f^ = 1), which satisfies the equation of parallel transport along the photon path. 

The electric field measured by a comoving observer is derived from the electromagnetic 
field tensor Ff"" = A^'f - A^'" as 

= = sR{[(a^'° - a°'^) + iw(aV - aV)]e*^} • (26) 

The length-scale over which the vector amplitude varies being much longer than the electromagnetic 
wavelength, the first term in the r.h.s. of Eq.(|2^) is negligible compared to the second. Thus, the 
direction of the electric field is indeed given by the direction of a*, in turn determined by /*. 

We shall calculate the projection of the electric field direction into the plane perpendicular 
to the x^-axis, with a = 1,2. Even though the photon paths deviate away from the polar 
axis in the presence of metric fluctuations, with typical deflections of order h, the projection into 
the plane perpendicular to unperturbed photon paths differs from the projection into the plane 
transverse to the actual photon trajectory by terms of order h? . 

The transverse components of the polarization vector verify the parallel-transport equation, 
which along a ray that propagates nearly parallel to the x^— axis towards the origin read: 

df"" 1 

-j^ = ^6''\hkc,Q-hbc,z + h^,^b-h^b,c)r . (27) 

Notice that, to the order considered, the difference between the affine parameter z and the 
x'^— coordinate is now irrelevant. 
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In the case of scalar static perturbations, /iafe,3 = h^^6ab is the only non-zero term in the r.h.s. 
of this equation. The solution to eq. (27) is such that 

r{z = 0) = [6\ + ^Ah\]f{zs) . (28) 

Since A/i^f, is diagonal, and are modified at the same rate, and thus there is no rotation of 
the plane of polarization with respect to the fixed coordinate grid. In other words: if 

X{z) = arctan (^JT^^ (29) 

and Ax = x{zs) - x(0), then 

Ax = (scalar) . (30) 



In the presence of a plane and monochromatic gravitational wave with wave-vector k = A;m, 
eq. (^) has instead non-diagonal terms 

^ = Z!^e*fc-(i-'^)[/i3'^(k)m, - hsb{k)m'' + (1 - ^^i)h\{k)]f» . (31) 
az 2 

Integration of this equation between the source and the observer leads to 

f-{z = 0) = [5% - u;e\ + ^Ah\]f\zs) (32) 

where, as in the previous section, A/i^fe = hab{zs,te) — hab{0,to) denotes the variation in the 
gravitational wave amplitude between the events of emission and observation of the photons, and 
a; = ^(1 -|- ^)A/ix coincides with the value for the image rotation induced by gravity waves. The 
angle that the polarization forms with respect to the -axis in this fixed coordinate grid thus 
changes by the amount 

Ax = + -(1 + y?)/S.hj^ sin 2(99 — x) + -^fJ-Ahx cos 2{ip — x) (tensor) . (33) 

The first term in the r.h.s, independent of the initial direction of polarization, coincides with the 
image rotation uj derived in the previous section. The other two terms, that arise from the presence 
of Ahab in eq. (|32|), coincide with the apparent rotation of a direction that subtends an angle x 
with respect to the x^-axis as a consequence of the local variation in the gravitational potentials 
between the source and observer locations (as we shall discuss further in the next section). 

Clearly, the rotation of the plane of polarization does not accumulate over distances longer 
than the wavelengths as light travels through a stochastic background of gravitational waves, 
since the effect is proportional to the change in the wave amplitude between the emission and 
observation events. The rotation of the plane of polarization produced by a stochastic background 
of gravitational waves of cosmological wavelengths and amplitude of order 10'^ is thus totally 
negligible. 



- 13 - 



4. Apparent rotations 

The rotation of the plane of polarization as light propagates across gravity waves derived in 
the previous section, or the lack of it in the case of scalar perturbations, were measured with 
respect to the fixed coordinate grid Notice, however, that due to weak lensing shear 

and rotation, be it of scalar or tensor nature, light originated in a point located at 9"'Zs, in a 
direction that forms an angle ^ = arctan(^^/^^) with respect to the x^-axis in the fixed coordinate 
grid, is seen by the observer at the origin as coming from a direction, in this same coordinate grid, 
at an angle ^ + = arctan(0^/^^), with 

A^ = w + 7sin2((/)-C) • (34) 

We recall that 7 is the shear intensity and (p its direction. A given fixed direction at the source is 
seen by the observer not just rotated by uj with respect to the absolute direction in this Euclidean 
grid, but is also affected by shear. In the case of scalar fluctuations, the apparent rotation of a 
fixed direction ^ induced by a single Fourier mode is 

= 7sin2((^ — ^) (scalar). (35) 

In the case of gravity waves, it is given by 

AC= u + l{l + fj.^)Ah+sm2{ip-C) + yAhy,cos2{ip-^) ^^^^ 
+7+ sin 2{ip — C) ~ 7x cos 2{ip — ^) (tensor). 

The first term in the r.h.s. of this equation is the rigid rotation uj. The next two terms are the 
effect upon a given fixed coordinate direction of the change in the gravitational potential between 
the emission and observation events. The rest is the additional apparent rotation due to the 
proper shear 7, given by eq. (p^). 

The question to address is what are the potentially observable consequences of apparent 
rotations of directions in the source relative to its plane of polarization. A possible approach is 
to consider sources with some structural, intrinsic direction. We can then answer the question 
of what happens to the apparent angle between the plane of polarization and that structural 
direction as light from the source propagates across metric perturbations. 

Consider sources with a significant degree of linear polarization and elongated shape, with 
approximately elliptical isophotes. For these objects we denote with ip the position angle on the 
sky (the orientation of the major axis of the image), and with x the direction of the plane of 
integrated linear polarization, after subtraction of the wavelength-dependent Faraday rotation 
along intervening magnetic fields. 

Weak lensing by metric perturbations produces focusing, shear, and rotation. Convergence 
amplifies sources without change in their shape, so it does not affect the observed position angle 
ip. Rotation, as described by to in eq. (12), accounts for a rigid rotation of the image, with the 
distance between any pair of points unchanged. Tensor perturbations rigidly rotate the position 
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angle ip by the amount oj. However, even when no rigid rotation occurs, as is the case for scalar 
perturbations, shear distorts a source image in such a way that an ellipse is still seen as an ellipse 
but with different ellipticity and orientation. 

Consider an elliptical source whose principal axis forms an angle ■0 with respect to the 
x^-axis in the source plane, = O'^Zg is taken to describe such an ellipse. We denote the length 
of the major and minor axes by o and b respectively. We define the ellipticity of the source 
as e = (o^ — 6^)/(a^ + 6^). When light from this object arrives to the observer distorted by 
convergence, shear, and rotation, the image parametrization j;" = {6°'i, + iIj°'i^)9^Zs in the observer 
sphere still describes an ellipse, but rotated an angle A-i/; and with a different ellipticity e', given 
by (if the distortion is small) : 

^ ^' l + 2(7/e)cos2((^- V) 



, _ V[e + 27 cos 2{(t) - V')]^ + 4[ea; + 7 sin 2((^ - ip) 



(38) 



1 + 2e7cos2(0 - V) 

Here (j) is, as in previous sections, the direction of the shear principal axis. If the ellipticity is 
larger than the shear, e ^ 7, we can approximate 

AV'«w + ^sin2(0-V) (39) 



(which of course coincides with eq. ( |34D in the limit e = 1), and the new ellipticity is 

e'«e + 27(l-e)cos2((?:)-'(/') ■ (40) 

In the opposite limit, of negligible eccentricity, e' ~ 27 and A0 « which simply means that the 
nearly circular source is distorted in the direction of the shear. 

We conclude that the angle between the plane of integrated linear polarization and the 
source axis, X ~ is observed to change, as light propagates across density inhomogeneities and 



gravitational waves, by the difference A(x — ip) as given by eqs. (30,B3|) and (37) 



Rotation and shear produced by a cosmological background of gravitational waves, at most 
of order 10"^, lead to a totally negligible apparent rotation between the plane of polarization and 
the source position angle. 

Scalar shear leads to an apparent rotation A(x — ^) = — A'i/', since A^ = for scalar 
perturbations. In sources with large ellipticity (e ~ 1), the apparent rotation is of the order of the 
shear itself. Such a small apparent rotation, at most of the order of a few degrees, is not likely 
to have significant observational consequences. Gravitational distortions due to weak lensing by 
large-scale structure will produce no observable effect on the polarization properties of extended 
radio jets in quasars. In the opposite limit, nearly circular sources appear distorted in the direction 
of the scalar shear, while their polarization plane remains in the same direction. Scalar shear 
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is thus able to mask any correlation between orientation and polarization that a nearly circular 
source might have. 

The most interesting regime may be the case of sources with ellipticities of the order of 
e « 0.1. These are the smallest ellipticities among radio sources that still present some significant 
(larger than 5%) degree of linear polarization. They can display an apparent rotation of the plane 
of polarization relative to their major axis of order 

A(x-V)«^ . (41) 

A scalar shear of 1% can change the observed angle between polarization and source orientation 
by 5 degrees. 



5. Effects of shear upon the correlation between source polarization and orientation 

The direction of the plane of polarization in radio galaxies and quasars is expected to trace the 
geometric configuration of the magnetic field in such objects, at least in sources with a significant 
degree of linear polarization. Discrepancies between the observed values and expectations based 
in models for the magnetic field in distant radio sources could eventually signal the presence of 
unusual effects additional to the standard Faraday rotation in the polarization properties of light 
as it propagates through the universe ( pirch 1982| , Carroll, Field fc Jackiw 1990| , Harari &; Sikivie 



1992, Nodland fc Ralston 1997, Carroll k Field 1997). 



It is possible that weak lensing by scalar perturbations explain the dispersion in the direction 
of integrated polarization of cosmologically distant radio sources away from the perpendicular to 
their major axis. We now discuss this possibility. 

Consider the data on 160 radio sources used in Carroll et. al (1990), with the corrections 
noted in Carroll and Field (1997). There are 89 sources with redshift z < 0.3 and 71 with redshift 
z > 0.3. The histogram of number of radio sources vs. X ~ V' ™ the case of nearby (z < 0.3) 
galaxies does not display a very strong correlation, although there is a narrow enhancement around 
X — tp ^ 0° and a broader peak at x — "0 ~ 90°. In the case of distant {z > 0.3) sources, the peak 
around x — "0 ~ 90° enhances dramatically, while there is no noticeabe peak around x ~ ~ 0°- 
It is likely that any correlation between polarization and source orientation be weaker in the case 
of sources only weakly polarized. Restriction of the sample to the case of sources with maximum 
polarization Pmax > 5% does indeed significantly enhance the two peaks in the histogram of nearby 
sources, while the peak around x — V-' ^ 90° ^oi distant sources becomes even more convincing. 

The most likely explanation for these observations is that the sample consists of two different 
populations. One population would consist of low luminosity sources, for which either there is no 
significant correlation between polarization and position angle, particularly in the case of weakly 
polarized sources, or have x ~ ~ 0°- The low luminosity population would be underepresented in 
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the high redshift observations, where the second class, of high luminosity sources with polarization 
almost perpendicular to their axis, would dominate ( |Carron &: Field 1997] , [Clarke et al. 198q ). 



We now put forward the following scenario, which appears to be consistent with present 
observations. Assume that most of the distant radio sources that have a significant degree of 
integrated linear polarization are polarized in a direction almost perpendicular to their major axis 
(say less than 5° away from it), while a small fraction does not posses a significant correlation 
between polarization and position angle. Scalar shear then acts stochastically upon the light 
emitted by these sources. The result should be an increase in the spread of the observed x ~ 
away from 90°, proportional to the rms value of the shear. The effect actually depends upon 
the source redshift and ellipticity. If the number of data points were sufficiently large one could 
search for a correlation between departures of x — away from 90° with the source redshift and 
ellipticity, and put to test this scenario. Present observations are not sufficient to make serious 
attempts in this direction, but we shall argue that they are at least consistent with the above 
scenario, if the scalar shear has rms value around 5%. Notice that this shear needs to be coherent 
only over small scales, of the order of the angular size of the source. 

Given that the number of data points is not sufficiently large to search for any correlation 
of the effect of shear with the source redshift and ellipticity, we shall estimate the effect to be 
expected in a population of radio sources with ellipticity equal to the mean in the sample, at a 
redshift equivalent to the mean redshift. Since shear effects do not significantly accumulate in 
the case of nearby sources, we choose to consider just the 51 sources with z > 0.3 and maximum 
integrated linear polarization Pmax > 5%. We plot the corresponding histogram in Figure 1. These 
sources have a mean value of x — V' equal to 90° with rms dispersion a = 33° , and mean redshift 
z = 0.85. Of these 51 sources, 41 have their angle of intrinsic polarization less than 45° away from 
the direction perpendicular to the source axis, while the remaining 10 appear to be more or less 
uniformly distributed at angles less than 45° away from the source axis, although the number is 
too small to draw significant conclusions. If we exclude from the sample these 10 sources, under 
the assumption that they may be predominantly members of a population with weak correlation 
between polarization and position angles, the remaining 41 sources have a mean value of x — ^ of 
91° with rms dispersion a = 19° and mean redshift z = 0.84. We have evaluated the ellipticites of 
all but 4 of them, for which we found the values of the major and minor axis quoted in Gregory at. 
al (1996). The mean ellipticity of the sample is e « 0.15. We put forward the hypothesis that the 
intrinsic dispersion oi x ~ 4' away from 90° is much smaller than the observed value, the observed 
dispersion being the consequence of the action of stochastic cosmic shear. 

We now estimate the amount of shear necessary to produce a mean square rotation a ~ 20° 
upon a population of sources with ellipticity e = 0.15, assumed to be located at equal redshift. 
In the limit 7/e << 1, where 7 is the rms value of the shear, assumed to act in all directions in 
the plane perpendicular to the photon paths with uniform probability, the results of the previous 
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section (eq. (|3^)) indicate that the rms value of the apparent rotation is 



1/2, 



~ / ' 



nl/2 



V2 e 



(42) 



A value = 0.35 (20°) corresponds to 7/e « 0.5, which implies a rms shear of 7.5% if the 
ellipticity of the sources is e = 0.15. Eq. (|39|), however, underestimates the apparent rotation when 
7/e is not significantly smaller than unity, as is frequently the case if the rms value is 7/e = 0.5. 
To better estimate the rms apparent rotation in a general case we should use eq. (|37D. Assuming 
a gaussian distribution of 7 around a zero mean, with rms dispersion 7, and uniform probability 
distribution in the shear direction 0, the rms rotation reads 

1/2 



2 /"Va 

~ / ' 

vr Jo 



2(7/e) sin 2(1) 



(43) 
In the 



,l + 2(7/e)cos2(/)y ^J^-n^ 

In the limit 7/e << 1, this expression coincides with the approximate value of eq. 
opposite limit, 7/e >> 1, it tends to 7r/(2\/3) (k. 52°), which corresponds to a uniform probability 
distribution for rotations between and 7r/2, as expected since in this limit the apparent rotation 
is given by "0 = i?^- 



The result of the numerical evaluation of eq. (43) is plotted in Figure 2. It indicates that 
= 0.35 (20°) when 7/e 0.3, which if e « 0.15 implies 

7 ^ 0.05 . (44) 

The statistical significance of the limited sample analysed, as well as the crude approach taken 
towards the redshift and ellipticity dependence of the effect, make this result only an approximate 
estimate. It seriously indicates, however, that the observed departures of the direction of intrinsic 
polarization away from the perpendicular to the source axis in distant radio sources are consistent 
with the hypothesis that their origin is the effect of shear with rms value around or below 5%. 
Much larger values of 7 would actually conflict with the observed correlation between source 
polarization and orientation. 

Other accurate measurements of the polarization properties of distant sources come from 



high resolution measurements of highly polarized local emission regions ( Leahy 1997 , War die et al. 



1997 ). In these observations, a very tight correlation is found between the plane of polarization 
and some structural direction (such as the direction of intensity gradients), in agreement with 
theoretical expectations. Shear by weak lensing effects would rotate the observed angle between 
polarization and that locally defined direction by the rms amount Aip « 7/\/2 only. 



6. Conclusions 



Measurement of the rotation of the intrinsic polarization with respect to the local orientation 
of each jet segment in extended radio jets has been used to probe the mass distribution in 
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intervening lensing galaxies ( [Kronberg et al. 1991 , Kronberg, Dyer, fc Roser 1996| ). In this paper 



we extended this method of gravitational "alignment-breaking" to the case where the change in 
the observed polarization properties of distant radio sources is due to the weak lensing effect 
produced by large scale density inhomogeneities. Indeed, we argued that cosmic shear can have 
a significant impact upon the observed angle x ~ '4' between the direction of linear polarization 
and some structural direction of radio sources. In extended radio sources (sources with large 
ellipticities e ~ 1) we found that weak lensing effects are negligible, since the change in the relative 
angle x~ ''P does not exceed the magnitud of the cosmic shear. On the other hand, radio sources 
with intermediate ellipticities, of the order of e ~ 0.1, can display a significant apparent rotation of 
the plane of polarization relative to their major axis due to gravitational distortions by large-scale 
structure. 

Our study provides an alternative approach to evidence the effects of cosmic shear coherent 
over small angular scales, of the order of the source angular diameter. This approach is 
complementary to methods currently used to detect cosmic shear coherent over the extension of 
galactic fields through the search of correlations in the ellipticities of the galaxies. 

Our main conclusion is that density inhomogeneities do not rotate the direction of linear 
polarization, but since they distort elliptical shapes, the angle between the direction of linear 
polarization and the source major axis observed from Earth differs from the angle at the source 
(after subtraction of Faraday rotation). In the limit 7 << e, where 7 is the value of the shear 
induced by the perturbations between the source and the observer, and e is the source ellipticity, 
the apparent rotation is of order A^/; ~ ^. This becomes very interesting in the case of sources 
with ellipticity of order e ~ 0.1, since it implies, in the presence of a shear of order 1 %, a rotation 
of around 5 degrees, already comparable or larger than the typical errors of a few degrees in 
Faraday rotation subtraction. 

We have shown that the observed correlation between polarization properties and source 
orientation in distant radio galaxies and quasars can be used to either constrain or eventually 
reveal the presence of cosmic shear. The number of present observations does not appear to be 
large enough to test with sufficient statistical significance whether the observed departures of 
the intrinsic polarization angle away from the direction perpendicular to the source axis may 
be ascribed to the effects of scalar shear. Present observations show that the largest fraction of 
distant radio sources have integrated linear polarization in a direction almost perpendicular to 
the source axis, with a dispersion around the mean of about 20°. We have argued that current 
observations are consistent with the hypothesis that this dispersion is caused by a shear with rms 
value around 5%, acting upon a population of radio sources that has a much smaller intrinsic 
dispersion. A rms shear much larger than 5% would induce rms rotations much larger than 20°, 
and would thus be inconsistent with the observed correlation between source polarization and 
orientation. With a larger number of observations available, the cosmic shear hypothesis could 
be tested through statistical correlations between the polarization vs. orientation angle and the 
source redshift and ellipticity. 
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Finally, weak lensing by a cosmological background of gravity waves is mostly of academic 
interest, since their effect does not accumulate with distance as that of scalar perturbations, 
and thus the expected distortions are completely negligible (of order 10~^). From an academic 
standpoint, however, it is interesting to emphasize the differences between tensor and scalar 
distortions. Gravity waves produce no focusing to linear order in the metric perturbation. 
They do induce shear plus rigid rotation of images, and also rotate the plane of polarization of 
electromagnetic radiation. 
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Fig. 1. — Histogram of number of sources vs. x~''P (polarization - position angle) for the 51 sources 
with redshift z > 0.3 and maximum integrated polarization Pmax > 5% 




Fig. 2. — Root mean square value of the apparent image rotation cr^ induced by a rms shear 7 
upon a source with elhpticity e, plotted as a function of 7/e 



